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The isotropic-to-nematic transition in an athermal solution of long rigid rods subject to a grav- 
£f~^ ■ itational (or centrifugal) field is theoretically considered in the Onsager approximation. The new 

' feature emerging in the presence of gravity is a concentration gradient which coupled with the 

f^**) | nematic ordering. For rodlike molecules this effect becomes noticeable at centrifugal acceleration 

£SJ , g ~ 10 3 -r 10 4 m/s 2 , while for biological rodlike objects, such as tobacco mosaic virus, TMV, the 

effect is important even for normal gravitational acceleration conditions. Rods are concentrated 
near the bottom of the vessel which sometimes leads to gravity induced nematic ordering. The 
concentration range corresponding to phase separation increases with increasing g. In the region of 
phase separation the local rod concentration, as well as the order parameter, follow a step function 
with height. 
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PACS numbers: 61.30.Cz, 64.70.Md, 61.25.He 
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Nematic ordering in a solution of long rigid rods has been studied theoretically in many papers, starting from 
the classical papers by Onsager and Flory 0. However, there is one aspect of this problem, which has never 
been considered, namely that this transition always occurs in a gravitational field. This field induces a concentration 
inhomogeneity within the volume where nematic ordering takes place. Such inhomogeneity should, in principle, change 
some of the characteristics of the liquid-crystalline transition. 

The dimensionless parameter associated with gravitational field is j3 = mgh/kT where m = mo — pv is the mass of 
Oh a rod corrected for buoyancy (jo and v are the density of pure solvent and the volume of the rod, respectively) , h is the 
height, g is the gravitational acceleration, T is the temperature and k is the Boltzmann constant. For h = 1 cm and 
I ' room temperature conditions, this gives the following criterion: the inhomogeneity due to the normal gravitational 
, acceleration (9.8 m/s 2 ) becomes important for molecular masses of rod more than 10 7 g/mol. Thus, for rigid rods 
7— I ■ made of common synthetic macromolecules (mo ~ 10 5 g/mol) this effect can be neglected. 

However, in at least two situations the effect of gravity on the problem of the liquid crystalline transition is important 
and experimentally relevant. First, for nematic ordering in solutions of high molecular weight rodlike biological objects 
(such as TMV or virus fd) 0,11 0,1113 m can be ver y lar § e ( m o = 4 x 1q7 g/ mo1 for TMV,0), and values of (3 > 1 
can easily be reached. Second, instead of normal gravitational acceleration, one can consider the acceleration in an 
ultracentrifuge which can be 10 4 -j- 10 5 times larger than ordinary gravity. For at least these cases, the investigation 
of the influence of gravitational field on the nematic ordering in the solution of rigid rods seems to be an important 
' problem. This problem is solved theoretically below in the Onsager approximation. 
• i— I ' 
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II. THEORETICAL MODEL 

The Onsager approach is based on a virial expansion of the free energy of the solution of rigid rods taking into 
account steric repulsion only. 

Let us choose the Cartesian coordinate system in such a way that the external field acts in the z direction and put 
the origin of coordinates at the bottom of the reservoir (z=0). Further, let us divide the volume of the vessel V into 
large number of identical layers aligned perpendicularly to the field in order that all particles in a given layer has the 
same gravitational potential. Here we use the following notation: dN z (0) = N z f z (SI) g?S1 is the number of rods in 
the layer z with axis directions within the small spatial angle dtt, N z is the total number of rods in the layer z, and 
£ = z/h denotes the dimensionless height. Let f z (SI) be one-particle orientational distribution function of rods in the 
layer z. The normalization for the function f z (SI) is written in the familiar form J f z (SI) dSl = 1. 

With such division of the volume of the vessel into very large number of layers, the local rod concentration C'(Q 
has the same value within the layer. Thus, for a given layer one can apply the traditional Onsager theory Q, M 
justified for homogeneous system. This theory is valid for dilute enough solutions of very long rods. 
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In this case, the local free energy of the layer labeled £ is written as 

™(C) 



kT 



lnC(C)+ / / (Jl,C) In (47r/(n,0)dn + C(C) / / (01,0/(^2,002 (7) dfiidfi 2 + //3C/(0,C)rfO (1) 



where the first term represents the entropy of a translational motion, the second term is the orientation entropy, the 
third term describes steric interaction of rods in the second virial approximation, and the last term is the average 
potential energy of a rod in an external gravitational field U ext (0 /kT = jjjrC = PC- 

To calculate the third term one assumes that the 7-dependent second virial coefficient a 2 is the half of excluded 
volume of two rods 0, thus, a 2 (7) = L 2 D |sin^|, where L and D are the length and diameter for long rigid rods and 
7 is the angle between directions of long axes. 

The free energy of the whole system is a sum of free energies of all layers. If the number of layers is large enough, 
the sum can be replaced by integral 

d(c' (C)^§-= f d(c' (0 (in d (0 + (7(0 + d (Op (0 + PC) (2) 



VkT 
where 



a(Q = J /(O,0M47T/(O,0)dO (3) 
is the orientational entropy of the layer £ and 

p'{C) =ff (Oi, / (o 2 , «2 (7) dfiidfia (4) 

is the second virial coefficient of interaction of two rods. 

To obtain the equilibrium distribution function we should take into account the possibility of the formation of a 
phase boundary between the nematic phase at the bottom of the vessel and isotropic phase on top. We denote the 
height of the boundary position in the vessel as £ = x, so the volumes occupied by the nematic and isotropic phases 
are V a — xV and Vi = (1 — x) V, respectively. With this, the free energy of the whole system eq. © becomes 

^(o^+/V;(o^#, (5) 



VkT J n " aw kT J„ * lw kT 

where 



^ = lnCl(0 + * a (0 + C a (0 Pa (C) + PC (6) 
is the local free energy of the nematic phase, and 

^ = in d i (0 + a. t (0 + d i (Opi (0 + PC (7) 

is the local free energy of the isotropic phase. 

To calculate the equilibrium distribution function / (Q, one should minimize the functional JSJ with respect to 
this function. The direct minimization of functional (JSJ) leads to a nonlinear integral equation, which can be solved 
only numerically 0, ^| . In the case where the volume of the vessel consists of two phases separated by a phase 
boundary, one should realize that / (f2, follows a step function with the variation of C with the function / a (0, i n 
the anisotropic part significantly different from that in the isotropic phase fi{C) = 1/4tt. To evaluate the distribution 
function f a (f2, m the nematic phase we apply an approximate variational method with a trial function depending 
on variational parameter a. 

Substituting this function in eqs. (J2J, ®, and minimizing with respect to a(C), we have the following equation 
for definition of variational function a(0 



(8) 
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However, the trial function proposed by Onsager [lj f(fl, a) — (a/ (Air sinh 0)) cosh (a cos O), where O is the angular 
deviation of a rod from the director, still leads to rather complicated integral equation. Therefore, following |!j for 
the sake of simplicity we used a trial function of simpler form 

tin a (0 j e- ;ii ^,0<e<^ ,» 

/, (fl, a(0) = ~rr { c,(Q(.-e)' _ _ ( 9 ) 



■in 



e - , f < 9 < ir 



with an approximate normalization (precise up to terms of order O (e ")). 

This trial function is suitable for approximate evaluation of <7 a (a(C)) in case of highly ordered state 



a a (a(C))= I / o (n,o(0)ln(47r/ o (n,a(C)))dnwlna(C)-l, (10) 
and the dimensionless second virial coefficient in the anisotropic phase 

Pa(C) = ts ^ i = - f a (n 1 ,a(C))f a (n 2 ,a(0)\ S m 7 \dQ 1 dn 2 ^^= (11) 



with notation = (^p (0/./. = ^ 2 -D (d sm 7l)i)j — jL 2 D; the value of 6 being equal to half of average excluded 

volume of two arbitrary oriented rods. 

The corresponding expressions in the isotropic phase are 

(Tj = and pi = 1. (12) 

In the above formulas the index i and a refers to the isotropic phase and the nematic phase respectively, and the 
angular brackets designate the average with respect to the isotropic distribution function f t (CI) = I/Att. 
Substituting the calculated values back in eq. (jHJ), yields the following expression for the function a(£): 



V a (C) = 



(13) 



where C a (£) = C a (£)6 is the dimensionless local rod concentration in the nematic phase. 
Thus, eqs. ®, and I© , © give 



roo(0 lnC;(C)+21n^^ + l + / 3C (14) 



fcT 
and 

^1=^^(0+^(0+^- (15) 
The chemical potentials in the phases can be also obtained 

Ha,i = V7 a ,i +Pa,i(() V a,i((), (16) 

where v a> i{C) — 1/C a i(0 i s the local specific volume and p o ,i(0 i s the pressure in the layer (: 

Pa,i(0 = 5 TTT = ( 6 a,i(0 a w ,^ (!') 

ft>o,i(C) v ' dC ai {C) 

The calculation of the pressure in the nematic and isotropic phases gives (compare with ref. 



p o (0 = 3(71(0, 



(18) 
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^(0 = ^(0(1 + ^(0). (19) 
Consequently, eqs. JEJ, JUJ, {T^J and dJ) give 

^ = In Cl (C) + 2 In + 4 + /3C (20) 

and 

^ = lnC.; (0+2CKO + 1 + /3C (21) 

In the equilibrium the chemical potential is independent on height and is the same in the both phases, thus one 
can obtain the equilibrium local concentrations in the phases: 

<(0 = 7^ T e-^ ) (22) 
h{x) 

where C a is the average concentration in the nematic phase, h(x) = J* e~^d£ is the normalization factor. 
Also 

Cl(0 = ^LW (2be-^ 1+ »>) , (23) 

where the function LW (x) corresponds to a solution of equation LW exp(LW) = x. For dilute solutions we can use 
the simple asymptotic form of this special function: LW (x) = x + 0(x 2 ). 
Thus, 

3(0 « y^e-^, (24) 

"ill: h{x) = J> N/O 

The equilibrium concentrations in the phases are determined by the following coexistence relations at the boundary 

Pa{x)=Pi{x), , 25) 

Substituting calculated values of chemical potentials (I2U|) . 121|) and pressures (|18|l . (|19|) with obtained concentrations 
in the phases (|22|l . (|23|l the coexistence relations (125[l is written as 

d 7TM e -7JM e l 1+ ; 7^ g J'- (26) 
In ^ + 2 In T + 3 = In ^ + 2^6"^, 

where C a = C a b and = C t b are dimensionless average concentrations in the nematic and isotropic phases, respec- 
tively. 

If the gravitational field is absent (g — > 0), eqs. (I26|) are reduced to corresponding equations for homogeneous 
system 
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III. OBTAINED RESULTS 

The numerical solution of eqs. (|26|l gives the values of average dimensionless concentrations of the nematic and 
isotropic phases, C a and Cj, coexisting at equilibrium. The phase diagram in the variables average rod concentration 
C - dimensionless parameter (3 is shown in Fig. ^ This diagram has three main regions. In the region labeled by 
letter /, entire solution of rods is isotropic (corresponding height of the boundary x = 0). In the region between 
the curves of coexistence Ci(/3) and C a (f3), the solution is separated in the isotropic and nematic phases with an 
interphase boundary between them (0 < x < 1). In the region designated by letter N, the entire volume of the vessel 
is occupied by the nematic phase (x = 1). 

This diagram reveals that gravity facilitates formation of the nematic phase (at least at the bottom of the vessel) 
and the region of phase separation becomes very broad even for rather low values of (3. 

Dependence of the local rod concentration on the height ( at fixed value of (3 is shown in Fig. The concentrations 
C a (C) and Ci(() obey the barometric distribution according to eqs. I|22|l and l|24|) . respectively. The concentrations of 
the nematic and isotropic phases in the boundary layer coincide with that in the absence of the field (C a (x) — 5.12 and 
C a (x) — 3.45 0)- This is the case because all the rods within a given layer exhibit the same gravitational potential. 
Thus, the rod concentration follows a step function with a jump at the phase boundary. 

The corresponding change with height of the order parameter S(Q) = J Pa (cos a(£l, £)) /(fi, C)d£l is shown in Fig. 
|3| The order parameter in the nematic phase is overstated due to the approximate trial function of the form 
(0. With increase of £ the order parameter decreases to the value corresponding to a nematic phase coexisting with 
isotropic one in the absence of gravity and then falls to zero. 

The Onsager approximation used in this paper is valid for low rod concentrations (volume fraction tp < 0.1). 
However, with increase of (3 the local rod concentration at the bottom of the vessel gradually increases. Thus, for 
high values of (3 barometric distribution (|22|l is no longer valid. 

To generalize the Onsager theory for the case of high rod concentrations one can use the Parsons approximation^^! 
which aims to improve the second virial coefficient by means of additional multiplier depending on mole fraction of 
rods, as well as some others generalizations (cf. refs. [la. HI . Il5t llq ) . Nevertheless, the calculations with nonanalitic 
distribution arising from such an approach are rather complicated and lead to additional integral equation. In most 
practical cases, except sedimentation in ultracentrifuge, the values of j3 are not too high (e.g. for TMV (3 is slightly 
above the unity), and traditional second virial approximation is quite justified. 

The position of the phase boundary vs. (3 for different values of total rod concentration is shown in Fig. These 
plots lead to the following conclusions. If the total rod concentration is low enough (i.e. the greater part of the vessel 
is occupied by the isotropic phase, C ~ Cj), the increasing gravity induces the isotropic-to- nematic transition and 
phase boundary shifts toward the top of the vessel. This process is observed until (3 ~ 1 and then the shift of the 
boundary stops, and the volume of nematic phase even slightly decreases (bottom phase is becoming denser under 
gravity; solid squares). If the total rod concentration is high enough (C ~ C a ) the nematic phase simply shrinks 
under gravity starting from the top of the vessel and the position of the phase boundary becomes gradually lower 
(solid circles). 

Furthermore, it is noteworthy to emphasize the important conclusion arising from the form of Fig. ^ The right 
branch of the plot rises very rapidly as (3 gets large, thus, remaining within the framework of general concepts dealing 
with spatially homogeneous phases one could suggest that the concentration of rods in the nematic phase should also 
rapidly increase with (3. However, the increase in the average concentration in the nematic phase is not as drastic as 
it follows from Fig. ^ The general reason is that the "rule of lever" cannot be applied for the present system, because 
we are dealing with spatially inhomogeneous phases. 

The average concentrations of nematic and isotropic phases, corresponding to phase separation, are shown in Fig. 
[5j This plot demonstrates that the average concentrations in both phases at fixed value of (3 do depend on the total 
concentration of rods and they do not coincide with the concentrations corresponding to the curves of Fig. ^ (dashed 
lines in Fig. [SJ. This is because (3 ~ h, where h is the total height of the vessel. Thus, the parameter (3 is different 
for separate phases and for the system as a whole. That is why the average concentrations of the phases lie within 
the region of phase separation shown in Fig. ^ 

IV. CONCLUSIONS 

Gravitational or centrifugal external fields facilitate liquid-crystalline transition at the bottom of the vessel and 
broadens the region of phase separation. This phenomenon should be noticeable for biological rod-like objects or 
common lyotropic molecules sedimenting in a centrifugal field. This seems to be an important problem which requires 
experimental investigation. 



Acknowledgments 

The authors thank Dr. S. Fraden who has drawn their attention to this unsolved problem. 



[1] Onsager L. // Ann. N. Y. Acad. Sci. 1949, V. 51, P. 627. 
[2] Flory P. J. // Proc. Roy. Soc. 1956, V. 234, P. 73. 

[3] Fraden S., Maret G., Caspar D. L. D, Meyer R. // Phys. Rev. Let., 1989. V. 63, P. 2068. 
[4] Fraden S., Maret G., Caspar D. L. D. // Phys. Rev. E, 1993, V. 48, P. 2816. 
[5] Nakamura H., Okano K. // Phys. Rev. Let., 1983, V. 50, P. 186. 
[6] Adams M., Fraden S. // Biophysical Journal, 1998, V. 74, P. 669. 
[7] Dogic Z., Fraden S. // Phys. Rev. Let., 1997, V. 78, P. 2417. 
[8] Starley J. P. // Mol. Cryst. Liq. Cryst, 1973, V. 22, P. 33. 
[9] Odijk T. // Macromolecules, 1986. V.19, P. 2313. 
[10] Kayser R. F., Raveche H. J. // Phys. Rev. A., 1978, V. 41, P. 53. 

[11] Lekkerkerker H. N. W., Coulon P., van der Haegen R., Deblieck R. // J. Chcm. Phys., 1984, V. 80, P. 3427. 

[12] Parsons J. D. // Phys. Rev. A., 1979, V. 19, P. 1225. 

[13] Lee S. D.// J. Chem. Phys., 1987, V. 87, P. 4972. 

[14] Hentschke R.// Macromolecules, 1990, V. 23, P. 1192. 

[15] Pre Du., Yang Y. // J. Chem. Phys., 1991, V. 94, P. 7466. 

[16] Chen J.// Macromolecules, 1993, V. 26, P. 3419. 



7 




Average Rod Concentration 



FIG. 1: Phase diagram in the variables (3 - C for athermal solution of long rigid rods. Here /3 = mgh/kT is the dimensionless 
parameter associated with external field, C is the dimensionless average rod concentration. Label N designates the nematic 
phase, I the isotropic phase while N + I corresponds to phase separation region. 




FIG. 2: Dependence of the dimensionsless local rod concentration C(C) on the dimensionless height ( at fixed value of /3. Open 
circles correspond to (3 = 1.1; solid circles correspond to (3 = 6.1; dimensionless total rod concentration is C = 3.4. 




FIG. 3: Dependence of the order parameter S on the dimensionless height £ at fixed values of (3. Open circles correspond to 
(3 = 1.1; solid circles correspond to (3 = 6.1; dimensionless total rod concentration is C = 3.4. 
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FIG. 4: The phase boundary position plotted as a function of j3 for different values of dimensionless total rod concentration 
C. The solid circles refer to C = 4.8; the open triangles to C = 4.1; the solid squares to C = 3.4. 
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FIG. 5: Average concentrations of rods in isotropic and nematic phases as a function of /3 for different values of total rod 
concentration C. Open circles correspond to C = 3.4, solid squares correspond to C = 4.8. 



